Early-time critical dynamics of lattices of coupled chaotic maps 

Philippe Marcq^ and Hugues Chate^ 

^ Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606, Japan 

^ CEA — Service de Physique de I'Etat Condense, Centre d' Etudes de Saclay, 91191 Gif-sur-Yvette, France 

(September 19, 1997) 

The early-time critical dynamics of continuous, Ising-like phase transitions is studied numeri- 
cally for two-dimensional lattices of coupled chaotic maps. Emphasis is laid on obtaining accurate 
estimates of the dynamic critical exponents 0' and z. The critical points of five different models are 
investigated, varying the mode of update, the coupling, and the local map. Our results suggest that 
the nature of update is a relevant parameter for dynamic universality classes of extended dynamical 
systems, generalizing results obtained previously for the static properties. They also indicate that 

J~>«. ' the universality observed for the static properties of Ising-like transitions of synchronously-updated 

^\ ' systems does not hold for their dynamic critical properties. 

0^■ 
o 

" ■ I. INTRODUCTION 

The last decade has seen considerable experimental, numerical, and analytical effort aimed at better understand- 
ing the sustained spatio-temporally chaotic regimes of large, homogeneous systems kept far-from-equilibrium by an 
^ \ external driving force. In particular, detailed investigation of a number of hydrodynamical flow regimes, including 
C^l ■ convection, shear flow, and crispation experiments, has led to a wealth of interesting insights into the properties of 
^^ ■ spatio-temporal chaos |l[ . Important pending questions concern the status of the asymptotic limit of long time and 
^— ^ , large system size, as well as the relationship which may exist between classical equilibrium statistical mechanics and 
^_j ' possible statistical descriptions of spatio-temporal chaos in this "thermodynamic" limit ||l|,|| . 

r — [ This article focuses on the critical behavior of models of spatio-temporal chaos close to second-order-like phase 

^\^ ■ transitions which occur in the thermodynamic limit. Theoretical work p] has suggested that phase transitions in 

"^^ [ generic non-equilibrium systems made up of locally interacting subunits belong to the universality class of Model 

^7' A B, for both static and dynamic critical exponents, provided that the order parameter is a non-conserved, scalar 

O quantity. Being based on coarse-grained Langevin descriptions, the approach developed in |g| overlooks the exact 

JL nature of microscopic time-evolution. Its conclusion also relies on the validity of assumptions generally associated 

jrt with the dynamic renormalization group formalism. 

'"^ Another significant contribution is that of Miller and Huse. In pj, they introduce a simple lattice dynamical system 

Ci ' with microscopic Ising symmetry (square lattice of locally-coupled, chaotic, odd maps), whose salient feature is the 

^ ' presence of a non-equilibrium continuous transition qualitatively similar to the ferromagnetic critical point of the 

two-dimensional Ising model. Ising-like transitions between spatio-temporally chaotic phases turn out to be a fairly 

common feature of coupled maps lattices (CMLs): they are observed for a variety of local maps, lattice geometries, 

and update rules l§,0j§!0l- However, contrary to the conjecture of and to early conclusions based on simulations 

of much smaller systems [B , careful analysis of finite-size data obtained from extensive numerical simulations of the 

Miller-Huse model shows that the corresponding phase transition does not belong to the Ising universality class [0,|| . 

The measured correlation-length exponent v = 0.89 ± 0.02 is significantly lower than vismg — 1, while exponent 

ratios (j/v and 7/1/ remain in good agreement with Ising values. Comparison with related models, in particular 

with transitions of sequentially-updated lattice dynamical systems, further indicates that synchronous update is the 

relevant parameter responsible for departure from Ising universality: keeping all other features of the Miller-Huse 

model unchanged, Ising static exponents are recovered as soon as sites are updated sequentially Isl. 

Our main objective is to extend previous work on static critical exponents to the dynamic critical properties of the 
Miller-Huse model. The dynamic critical exponent z, which quantifies the algebraic divergence of coherence times at 
criticality, is known to be sensitive to parameters otherwise irrelevant for static exponents, such as the existence or 
absence of macroscopic quantities conserved under time-evolution B . In addition, the nature of update is a relevant 
parameter for dynamical universality classes of Ising systems: synchronous update of clusters of spins yields distinct, 
significantly lower values of the dynamic exponent z than is observed for standard sequential or checkerboard update 
||10| . Ignoring both the conjecture of [^ and the numerical results of g], one may thus naively expect phase transitions 
of synchronously and sequentially updated CMLs to be characterized by different dynamical exponents. Here, we want 
to confirm, for the dynamic properties of Ising-like transitions of lattices of coupled chaotic maps, the relevance of the 



^ 



mode of update already discovered in [|| for their static properties. Similarly, the static universality class observed 
for synchronously-updated models is revisited from the point of view of their dynamical properties. 

Compared to Ising systems, the measurement of static critical exponents turns out to be significantly more resource- 
consuming in the case of CMLs, due in particular to the presence of unusually large corrections to scaling ||^ . Moreover, 
extracting reliable values of dynamic critical exponents from direct simulations is a notably difficult task, even in a 
priori simpler cases. Despite much numerical effort, the value of the dynamical exponent z of Model A in dimension 



d = 2 remains somewhat controversial (see ]ll[ | for a survey of work done prior to 1993, and |12| for a more recent 
review). The methodology we apply here to phase transitions of CMLs is based on recent theoretical work by Janssen 
et al., which proves the existence of a new universal regime in the early critical dynamics of systems starting from 
non-equilibrium (e.g. completely disordered) initial conditions [^^|jlj]. This regime, termed "initial critical slip" or 
"universal short-time behavior" in the literature, is characterized by a new non-trivial exponent 9' , unrelated to the 
usual static and dynamical exponents. The dynamical exponents 9' and z can be readily obtained from the initial 



scaling properties of observables of finite-size systems, as shown analytically in |15|, and first implemented numerically 
in ||l^. Unlike standard methods, this procedure is nearly free from the difficulties associated with critical slowing 
down at the transition point: useful simulation times (T r^ 10^ — 10^) are typically much shorter than the finite- 
size coherence time-scale t^ ~ L^ ■ Statistical accuracy is ensured by ensemble-averaging over a large number of 
independent realizations. Thanks to high numerical efficien cy, goo d agreement on the value of critical quantities such 
as 9' and d/z — 9' has been already reached for Model A | [l7| , p^ , p9[ . This makes comparison with other systems easier, 
and opens the way to an investigation of the relative universality of z and 9' which, based on the theoretical work of 
Janssen et al |l^Jl^, are expected to depend on the same relevant parameters. 

This article is organized as follows: current understanding of early-time critical dynamics is briefly reviewed in 
Sec. pi The methodology we follow closely parallels that used by Okano et al. for of the two-dimensional Ising model 
with heat-bath and Metropolis algorithm |l9]. The same procedure is used throughout, thus allowing meaningful 
comparison between exponents obtained for different CMLs, as well as with exponents of Model A. First, the dynamic 
critical properties of the Miller-Huse model, a lattice dynamical sys tem w ith synchronous update, are investigated in 



Sec. in. The model and its phenomenology are introduc ed in Sec 



the critical exponents 6' and z are next described in Sec. HI B and 



III A| . Simulations pertaining to the measure of 
III C respectively. In Sec. |^, we investigate the 



role played by the type of update for the dynamic critical properties of Ising-like transitions, in order to extend its 



relevance, already established in [B| at the static level. In Sec. IV A, we first consider a sequentially updated model 



introduced in |q], which, according to previous numerical results, belongs to the Ising universality class for static 



critical exponents. In Sec. IV B, we turn to Sakaguchi's model [B|, a CML with checkerboard update whose static 
critical exponents are known exactly to be equal to those of the Ising model. Next, we consider, in Sec.|y|, various 
synchronously- updated models to investigate whether the universality of the static critical properties of their Ising-like 



transitions extend to their dynamic exponents. Our results are summed up and discussed in Sec. VI 



II. EARLY-TIME CRITICAL DYNAMICS AT SECOND-ORDER TRANSITIONS 

In order to measure the dynamic critical exponent z from numerical simulations of finite-size systems, most methods 
considered until a few years ago made use of the so-called nonlinear relaxation regime, by, e.g., looking at the decay 
of the system's time-dependent magnetization M{t) according to: 

M{t)^t-'^/^'''\ (1) 

or similar relations involving higher-order moments. This regime was generally believed to be relevant within the time 
interval 1 <S^ t <S^ II, whereas finite-size linear relaxation eventually prevails beyond t^ = L^ , where the magnetization 
decays exponentially: M{t) ~ exp(— t/t^)- 

A point overlooked until the work of Janssen et al. |]l3| is the importance of initial conditions. Suppose that we 
start from disordered non-equilibrium initial conditions (magnetization is zero or very close to zero) with very short 
initial correlation length, and quench the system to its critical point. One qualitatively expects fluctuations to be 
negligible at first: the system is then mean- field-like. Since the mean-field critical temperature is usually larger than 
the actual critical temperature, the system is in its ordered phase, and the magnetization (and correlation length) 
will want to grow. This accounts for initial magnetization growth. There is of course a crossover point, after which 
the system's behavior reverts to the usual (relaxational) behavior of Eq. |l|. 

This qualitative idea has been formalized, and the influence of initial conditions on renormalization-group trans- 
formations investigated in detail for bulk systems |0,n. Let m^ be the initial magnetization at time t — Q. This 
field gives way to a new scaling index xq independent from already known ones (both static and dynamical), and to a 
timescale io within which a new universal scaling regime sets in. The new exponent 9' is universal in the same sense 



as the ususal dynamic critical exponent z, since it was obtained within the same formaUsm. For t,„ic < t < to, and 
mo small enough, the magnetization grows as a power law: 

Mit) ^ moi'', (2) 

where 9' = {xq — fijv)! z. The microscopic time imic is the time after which macroscopically correlated regions form, 
i.e. regions large compared to the microscopic length scale, in this case the lattice constant. The time-evolution of 
observables for t < tmic is nonuniversal, and depends on microscopic features of the model. The crossover time to is 
obtained by matching Eqns. (^ and (g): 

to ^m^^° \ (3) 

and diverges in the limit of zero initial magnetization toq. In that case, the nonlinear relaxation regime is not observed 
in the bulk. 

Then, finite-size scaling theory was introduced by [ p5[ . We will need it for interpretation of numerical experiments. 
The scale invariant expression reads, for a system of finite size L, and the fc**^ moment of the order parameter, at the 
critical point: 

Af W {t, L, mo) = h-^^'^'M^^'^ {h-H, b-^L, 6^°mo) , (4) 

where & is a scaling factor and M^*^' is a universal function, independent of microscopic details of the system. Important 
point: the exponent z in Eq. ^ is the same as the usual one Q. Choosing the arbitrary prefactor equal to 6 ^ t^^^, 
one obtains: 

M W (i, L, mo) =*-'=''/'' AfW(tAL,t/to) , (5) 

where to = ttiq ^° and tL — L^. For a finite-size system and evolution times t < to,tL, one obtains: 

Af(i)-moi^', (6) 

for small values of mo [Q . This allows to measure 9' directly. Then, assuming that the value of P/v is already known, 
z can be obtained thanks to the relation |13] : 

M(2) (t) - t'^ with C = (^ - W/v)/z , (7) 

as used for Model A in |1^ . Finally, careful renormalization group analysis leads to the following scaling form for the 
order-parameter time correlation function A{t) = {m(t)m{0)) (cf. explicit derivation in [Q) 

A{t) - t-^ with 6 = d/z-9' , (8) 

where the space dimension is denoted d. 

Note that finite-size scaling relations may also be used in order to measure z and P/v p(]|] : 

U{t,L) = U{b%bL), 

M(^'){t,L) = b^'^/''M^^Hb%bL). ^^' 

We choose not to, since finite-size effects seem to be either negligible, or easily controlled in cases relevant here (see 
below). 

To conclude this section, we briefly review recent work on the critical dynamics of Model A. The relevant exponent 
values are gathered in Table I. The exponent 9' has been measured twice according to Eqn. ^ first by Grassberger |l|] 
{9' — 0.191(3), heat-bath dynamics), then by Okano et al. |19] {9' — 0.194(4), heat-bath and Metropolis algorithms), 
thanks to slightly different methods. We use a conservative combination of the two estimates as our reference value: 



^ModelA 



0.193(5). (10) 



Excellent agreement has also been reached for the combination S = d/ z — 9' , obtained from Eq. g, between the early 
measures of Huse and of Humayun and Bray pTJ I {S = 0.74(1), heat-bath algorithm) and a recent confirmation by 
Okano et al. [|9| (S = 0.739(5)). Our conservative estimate is thus: 

^ModclA = 0.74(1). (11) 



The case of the dynamic critical exponent z is more delicate. Estimates using methods derived from the theory of 
early-time critical dynamics vary between 2.155(3) (heat-bath) ||l£], 2.137(11) (Metropolis) [g9|, and 2.143(5) (heat- 
bath) Ig^l- These estimates are somewhat lower than the currently accepted value z = 2.165(15) ||l^, obtained from 
both series expansions (z = 2.165(15), data from [2^ reanalized by Adler, see [Q) and from a number of direct 
simulations of very large systems: 2.165(10) (nonlinear relaxation ^M)^ 2.172(6) (damage spreading flq]), 2.160(5) 
(non-linear relaxation ||22|] ). Since the latter generally correspond to significantly better statistics and larger system 
sizes, we choose 

^ModelA = 2.165(15) (12) 

as our reference value. It leads to the combination C, = {d — 2(3/v)l z — 0.808(6), for (3/v = 1/8, in reasonable 
agreement, within error bars, with the value obtained from Eq. in |19|: C = 0.817(7). 



III. DYNAMIC CRITICAL EXPONENTS OF THE MILLER-HUSE MODEL 

A. The model 

Recently, Miller and Huse introduced a CML designed to be a simple non-equilibrium Ising-like model p| . Its local 
map, which provides the "reaction" part of this reaction-diffusion lattice dynamical system, is an odd, piecewise-linear, 
chaotic map of the real interval [—1,1]: 




(13) 
-3a; + 2 if \ <x< 1. 

The constant absolute value of the slope being equal to three, its Lyapunov exponent is positive and equal to In 3. 
For the simple case of a two-dimensional square lattice, the evolution rule reads: 

xiy = (1 - 4g) fixl) + g (/(x^i,,) + /(4,-i) + fi^Uu,) + /«,+i)) > (W) 

where t denotes the (discrete) time, and the subscripts the position on the lattice. The nearest-neighbor coupling 
constant g can vary between and 1/4. In the following, all numerical calculations are performed on square arrays of 
linear size L with periodic boundary conditions. 

Since f{x) is an odd function of x, discrete spin variables can be defined in a natural fashion: 

<^. =sign(4^.)e{-l,l}. (15) 

Next, the fluctuating, space- averaged magnetization is defined as: 

1 
l2 



< = l^E<r (16) 



ij 



In fact, one can also use a definition of the "magnetization" based on the original continuous variables xl .-. This does 
not alter significantly the statistical results, as we mention in the following. 

Increasing the coupling constant g, the only control parameter in this system, an Ising-like phase transition takes 
place from a disordered phase with zero average magnetization at weak coupling to an ordered phase at strong coupling, 
where the spins tend to be aligned with each other. The order parameter is the magnetization Ml = (I"^lI)i where 
the brackets represent in practice (long) time- averages (ergodicity is assumed). 

Note that chaos is extensive in this system B, and that dynamical quantifiers, such as the Kolmogorov-Sinai 
entropy, seem to be insensitive to the onset of long range order at least for the finite size systems for which these 
calculations can be made. Only one lengthscale, the correlation length ^, diverges in the thermodynamic limit. This 
proves that the transition exists in the thermodynamic limit, as corroborated by the applicability of finite-size scaling 
arguments. Such an analysis allows to measure the standard static critical exponents, and, in particular, the deviation 
of the correlation-length exponent v — 0.89(2), from the Ising value i/ = 1. g] 



B. Measure of 6' 

Here we want to look at the short-time dynamics of carefuhy-prepared initial configurations with a given initial 
magnetization mo. They are generated easily by the following procedure. For mo = 0, assign real random numbers 
{xij) uniformly distributed on [0,1] to L^/2 randomly chosen sites of the lattice. Assign then the opposite values 
(— Xij) randomly to the remaining L^/2 sites: the total magnetization is then exactly zero. In order to obtain a 
small, but non-zero initial magnetization, implement the same procedure, but based on {L'^ — K) 12 randomly chosen 
sites. Then set the value of the K other sites to, e.g., x = 1. The magnetization is thus equal to mo = 2K/L^ . We 
have checked that this particular choice does not influence the scaling properties described in the following. Choosing 
X = 1 possesses the advantage that the initial magnetization has the same value whether considering discrete spins 
a\ or the original continuous variables x\ ■ . 






(a) 




FIG. 1. Snapshots of the coarsening process: positive (resp. negative) spins are represented by a black (resp. white) pixel on 
a two-dimensional grid. This simulation of the Miller-Huse model is performed at the critical point Qc = 0.20534, for an initial 
magnetization mo — 2.4 10~^ {mo — 2K/L^,K = 2,L = 128). Snapshots (a), (6) and (c) correspond to evolution times t — Q 
(completely disordered initial condition), t — t^ic = 5 (build-up of macroscopic correlations) and t — T — 128 ("initial critical 
slip" regime). 

The value of the critical coupling strength gc was previously obtained according to Binder's method. We use 
gc — 0.20534(2) Q. As predicted in jl^], a regime of initial growth of the magnetization is observed, as well as 
the crossover toward nonlinear relaxation for large enough initial magnetization m-Q. The corresponding coarsening 
process is illustrated in Fig. |l|. For measurement purposes, we use K — 2,4, 6, 8, 10, for sizes ranging between L — 16 
and L = 128. The duration of a run is T = 128. In such conditions, no crossover is observed to the nonlinear 
relaxation regime, since t^ ^ T. Thanks to the good quality of our data, the value of the microscopic time imic = 5 
can be obtained by simple visual inspection (Fig. 0). This relatively small value is comparable to what has been 
observed for the two-dimensional Ising model |19| . Scaling of the magnetization versus time is observed over the time 
interval 5 < t < 128. This corresponds to the initial critical slip regime. 

The exponent 9' is measured thanks to a linear fit in log-log scale over the interval imic 1^ t ^ T. We checked that 
using larger values of imic and/or T does not alter the estimate. Note also that using values of the critical coupling 
outside the confidence interval gc = 0.20534(2) does not lead to an improved quality of fits: this confirms the validity 
of estimates of the critical coupling obtained in ||^ . 

Ensemble averages are performed over 512000 realizations for L < 64, 128000 realizations for L — 128. Statistical 
errors are estimated by comparing the exponent values obtained for five different initial magnetizations mo = 2K/L^, 
K — 2,4, 6, 8, 10. Error bars take into account the uncertainty on gc- Note that the corresponding values of mo are 
much smaller than those used by Okano et al., who needed to extrapolate exponent values obtained for small, but 
finite initial magnetization to the limit niQ — 0. Our procedure is similar to that used by Grassberger [p8|, since no 
extrapolation is needed. The exponent values thus measured for L = 32, 64, 128 are respectively 9' = 0.148(2), 0.142(4) 
and 0.148(7). Statistically equivalent values are obtained when considering the magnetization based on the continuous 
variables. Within error bars, no finite-size effect is observed for L > 32. Our global (conservative) estimate is: 



9}^H = 0.146(9). 



(17) 



Note that this result is not consistent with the accepted value for the critical dynamics of Model A: 6' = 0.193(6) 
obtained by similar methods and with a similar statistical quality in |1§| , |19|] . 
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FIG. 2. Measure of the dynamic critical exponent 9' for the Miller-Huse model. We plot in log-log scale the magnetization 
M{t) vs. time t measured at the critical point g ~ gc = 0.20534. Clear scaling is observed in all cases for t > tmic ~ 5. The 
solid lines in both graphs correspond to a slope equal to 9' — 0.146. Graph (a): the system size L = 128 is fixed, the initial 
magnetization mo = 2K/L^ varies between mo = 2.4 10~* and 1.2 10~^, for 5 distinct values oi K = 2, 4, 6, 8, 10 (from bottom 
to top). Good agreement between slopes corresponding to different values of K suggests that the limit mo -^ is satisfactorily 
approximated. Graph (6): the initial condition K = & is fixed, for system sizes L — 16,32,64,128 (from top to bottom). 
Finite-size effects are negligible. This suggests that the infinite-size limit is satisfactorily approximated. 



C. Measure of z 



A first way to measure exponent z is given by Eq. [7|. The experimental conditions are similar to those mentioned 
in the last section, but for an initial magnetization equal to zero (toq = 0, i^ = 0). Simulations were performed for 
L = 16, 32, 64, 128, T — 128, ensemble averages performed over 512000 realizations for L < 64 and 128000 realizations 
for L = 128. Log-log plots of the second moment M'-^'> vs time show very good scaling, and lead to estimates of the 
combination of exponents ( = (d — 2j3/v)/z = 7/(4z), assuming that [3/v = 1/8 for d = 2, as imphed by Q. In order 
to avoid interferences with the current experimental uncertainty on (i/v — estimated to be f3/v = 0.131(6) in [gj — , 
we will work with C, and convert into z as late as possible. Longer runs (up to T = 1024) were performed for large 
system sizes (L = 128), with lesser statistical accuracy. This allowed to check that the exponents measured do indeed 
correspond to the asymptotic regime (cf. Fig. |^). 

Here, determining the microscopic time imic requires additional effort, when compared to the previous case. We 
use a method introduced by Okano ct al. |jl9]. Local exponents C(i) are first measured from "local" fits limited to an 
interval of time [t, t ~\- iioc]. The microscopic time imic is defined as the time beyond which C,[t) becomes stationary, 
within statistical fluctuations. We find fmic ~ 30, for tioc = 15 and all sizes L considered. This value does not depend 
on the value of tioc for tjoc large enough. The microscopic time measured here is much larger than the one estimated 
for the scaling of M{t) at small vtiq, as had already been observed for the Ising model [|l9| . 

Asymptotic values of Q are then obtained, as before, from a global linear fit performed in log-log scale over the 
interval tmic 1^ t < T — 128. As before, we checked that values obtained do not change, within error bars, for 
15 < tmic ^ 50. Again, these values were independent on whether the discrete spins or the continuous variables were 
used to calculate the magnetization. 

In this case, finite-size effects are sizeable and can be well controlled: we obtain values of C(-^) for sizes 16 < L < 128 
which are monotonously decreasing and seem to converge. Statistical errors are evaluated from a comparison of 
exponent values measured for three coupling strength values in the confidence interval [0.20532,0.20536]. As before, 
the quality of fits does not increase when varying gc outside this interval. In order to evaluate the rate of convergence, 
we use the Ansatz: 



C(L) - C(oo) ^ L- 



(18) 



where C(oo) is the desired infinite-size exponent. Although algebraic relaxation toward an asymptotic value seems 
natural in the context of critical phenomena, we are not aware of any theoretical justification for Eq. 18. The validity 
of this phenomenological Ansatz is confirmed by our data (Fig. 0), which is not compatible with, e.g., exponential 
relaxation. Optimizing linear fits in log-log scale of C,{L) — C,{oo) vs. L yields the following estimate: C,{oo) — 0.839(3) 
(and, incidentally, w = 1.4). The Ising value 0.817 is not compatible with our data plus Ansatz ^^ (see inserts of 
Fig. g). Note that we use here the numerical results of [|9| as a reference value for Model A. The discrepancy between 
our data for the Miller-Huse model and the currently accepted value |1^ for Model A (^ = 0.808(6), see discussion in 
Sec. ||) is even larger. 

Using the theoretical, exact Ising value for j3/i> [g|, our estimate yields z = 2.082(4), at variance with the Model A 
exponent, irrespective of the method used to estimate it, whether it is initial critical slip |T9,EQ] or standard methods 
||l2| . Using the measured value (i/v = 0.131(6) ||], one obtains: 



^MH = 2.07(2), 



(19) 



a conservative estimate which we endorse. 
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FIG. 3. Measure of the dynamic critical exponent z for the Miller-Huse model. Graph (a): we plot in log- log scale the 
second moment M^^' (t) vs. time t at the critical point g = Qc = 0.20534, for an initial magnetization mo — 0. System sizes are 
L = 16, 32, 64, 128 (from top to bottom). Scaling is observed for t > tmic ~ 30. The solid line corresponds to a slope equal to 
(^ — {d — 2j3/u)/ z = 0.839. The curves corresponding to L = 64 and L — 128 (T — 128) are undistinguishable. This suggests 
that finite-size effects become small for L > 64. Data obtained for L = 128, T = 1024, 32000 realizations (triangular symbols) 
shows no evidence of cross-over to a distinct time-asymptotic regime. Insert: we plot the ratio M"'(f)/f'' vs. time t in log-lin 
scale for the values ( — 0.817 (Model A, top curve) and ( — 0.839 (best fit for extrapolated value in the infinite-size limit (see 
(b)), bottom curve). The system size is L = 128, T — 128. This graph shows that our (finite-size) data is not compatible with 
exponent values expected for Model A. Graph (6): we plot the finite-size exponent C^(L) vs. system size L, as obtained from 
linear fits of the data presented in Graph (a) for 30 < i < 128. The infinite-size estimate ( = C^(oo) = 0.839 (dashed line) is 
derived from the relation C(I/) — 0.839 + 3.9L^^'^ (solid line). The insert shows a log-log plot of (,{L) — C(oo) vs. system size 
L for the two values ( — 0.817 (Model A) and (" — 0.839 (measured). The latter value optimizes the quality of a linear fit of 
log(C(^) ~ vs. log(L). This plot indicates that our data is not compatible with exponent values expected for Model A when 
extrapolated to the infinite-size limit. 



Confirmation of the previously measured values may in principle be obtained from the scaling behavior of the 
correlation function (Eq. p[), with an exponent 5 = d/z — 9' . Previously measured values lead to d/z — 6' = 0.81(1). 
Even though the runs used are the same as for the second moment M*^^\ in practice, scaling is not satisfactory. 
In particular, local exponents Siocit) do not converge to stationary values. Very strong finite-size corrections to the 
dominant scaling are present, which preclude any effective measurement of 6. However, the observed behavior is 
compatible with large-size, long-time convergence to the above value 5 = d/z — 6' = 0.81(1). Direct, reliable estimates 
for S remain beyond our numerical resources. 



In conclusion, the above study of the Miller-Huse model first shows that the regime of initial critical slip exists also 
for (some) deterministic systems, with the same phenomenology as for spin systems. Quantitatively, our simulations 
lead to the conclusion that the Miller-Huse model does not belong to the dynamic universality class of Model A. 
Since the mode of update is the relevant parameter explaining departure from the Ising universality for static critical 
exponents in this model [p|, one would naturally like to know whether or not this is also true for dynamic exponents. 
Section IV deals with this question, with the study of two models with respectively sequential and checkerboard 
update. 



IV. UPDATE RULES 



A. Sequential update 



The model considered here is identical to the Miller-Huse model except for one point: the update rule. Evolution 
rule E4 is replaced by: 



^1 = (1 - 45) /«,-) + g {lixth) + /(^■+ii + K4+1,,) + /«,+i)) • 



(20) 



Sites are updated one at a time, in sequential order from the top, leftmost site to the bottom, rightmost one, as in: 



(1,1) 
(1,2) 



(2,1)^(3,1) 
(2,2)^(3,2) 



(1,L)^(2,L)^(3,L) 



(i,l)- 
(i,2)- 



(21) 



where arrows indicate the order of update between sites of indices (i, j). Consequently, boundary conditions are helical 
in this case. (Note that this mode of update was termed "asynchronous" in [g|.) A continuous transition similar to the 
Ising ferromagnetic point occurs in this system too, albeit for significantly lower coupling strength gc = 0.11255(5), 
as measured in B thanks to Binder's method. Its (measured) static critical exponents are compatible with the static 
Ising universality class P/u = 0.117(12), v = 1.02(7) §]. 
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FIG. 4. Measure of the dynamic critical exponent 9' in the case of sequential update. The magnetization M{t), measured 
at the critical point g = Qc = 0.11255, is plotted vs. time t in log-log scale. Scaling is observed in all cases for t > tmic = 5. 
The solid line in both graphs corresponds to a slope equal to 6' — 0.06. Graph (a): the system size L — 128 is fixed, the initial 
magnetization mo = 2K/L^ varies between mo = 2.4 10~^ and 1.2 10""^, for values oi K = 2,4,6,8, 10 (from bottom to top). 
Graph (6): the initial condition A' = 6 is fixed, for system sizes L — 16, 32, 64, 128 (from top to bottom). Finite-size effects are 
negligible. 



The initial conditions to measure the early-time critical properties of this transition are prepared as for the Miller- 
Huse model. The observed phenomenology is the same, for a comparable microscopic time imic '^ 5, as estimated 
visually (Fig. ||). The value of imic is confirmed by the method described in Sec. [IIC| (local exponent). Finite-size 
effects are negligible for L > 32. Measurement of the dynamic critical exponent 9' is based on sizes L = 32, 64, 128. 
The number of realizations over which ensemble-averaging is performed is the same as before, 512000 for L = 32, 64, 
and 128000 for L = 128 (see Fig. |). Our final estimate is 

^Sequential = 0.06(2). (22) 

Error bars take into account both statistical errors and the uncertainty due to the error bars on the location of the 
critical coupling strength. This is clearly not compatible with either the value obtained for the Miller-Huse model 
{9' = 0.146(9)) or with the Model A value {9' = 0.193(5)). At this stage we have no particular understanding of 
the low value of 9' in this case. Finally, note that our data is equally well fitted by a logarithmic time-dependence: 
M{t) ~ log(t), for which no theoretical justification is at present available. 
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FIG. 5. Measure of the dynamic critical exponent z in the case of sequential update. The second moment of the magnetization 
-^^ (t) (graph (a)) and the temporal autocorrelation function A{t) (graph (6)) are plotted in log-log scale vs. time t, as measured 
at the critical point g — 0.11255, for zero initial magnetization mo — 0, up to time T — 2048. The system sizes considered are 
L = 16,32,64,128. All curves superpose for M^'^'{t), and the last three for A{t). The microscopic time tmic can be roughly 
evaluated to tmic = 0(1000), especially from the inserts, where we plot in log-lin scale the local exponents Cioc(i) and (5ioc(i) 
vs. time t, computed for L — 64, over time intervals of duration iioc ~ 100. Graph (a): the (^ values corresponding to Model 
A (^ = 0.817, solid lines) and the Miller-Huse model (^ = 0.839, dashed lines) are indicated. Graph (b): the bottom and 
top curves respectively correspond to sizes L — 16 and L = 32,64, 128 (superposed). The 5 values corresponding to Model A 
{5 — 0.74, solid lines) and the Miller-Huse model (^ — 0.81, dashed lines) are indicated. Note that fluctuations of Cioc(i) and 
Sioc{t) observed for t > imic encompass the values for both models. 



Measuring the exponent z is much more difficult, since the microscopic time (or rather the time beyond which 
corrections to dominant scaling vanish and the asymptotic regime sets in) is very large, of the order of iniic ~ 1000. 
This is clearly shown in Fig.|5|, and in particular from the plots of the local exponents vs. time in the inserts. 
Ensemble averages are computed over typically 64000 independent runs, for a simulation time of T = 2048. Estimates 
of exponents must be derived from intervals much smaller than one decade (typically 1000 < t < 2048). Linear fitting 
in log-log scale is thus impractical. For that reason, we assume that the plateau observed for t > t^ic in plots of 
the local exponents correspond to the asymptotic value, and evaluate error bars from the variation of local exponent 
within that interval (FigJ^). We obtain: 

Cscqucntial — 0.83(6), 
^Sequential = 0.78(8). 



(23) 



The estimated value of C, close to that of the Miller-Huse model, leads to z = 2.12(15) (assuming fS/v ~ 1/8), with 
error bars large enough to include the uncertainty on the exponent of both Model A and Miller-Huse model. The 
estimated value of 6 is closer to that of Model A, but is too imprecise to be exploited. 



We are thus unable to give an estimate of z which would decide between the values for the Model A, the Miller- 
Huse model, or an eventual third number. It is clear, though, that our estimate of exponent B' {9' = 0.06(2)) is not 
compatible with the values of either Model A {6' = 0.193(5)) or the Miller-Huse model {9' = 0.146(9)). 

The above results confirm, at the dynamic level, those obtained in [|| for the static critical properties of Ising- 
like phase transitions in systems made of coupled chaotic maps: the mode of update is relevant. In |p|, the static 
exponents of the sequential-update model studied above were measured to be, within numerical accuracy, those of 
the Ising model. The above results show that these two models have different dynamic exponents (at least 9'). Is 
this due to their mode of update and/or to other factors? In the next subsection, we introduce Sakaguchi's model, a 
system of coupled chaotic maps designed to have exactly the static exponents of the Ising model, in order to further 
assess the relevance of the mode of update and that of the nature of the model for the dynamic scaling properties of 
Ising-like transitions. 

B. Checkerboard update 

A few years ago, Sakaguchi introduced a system of coupled Bernoulli maps with an exponential coupling scheme 
and checkerboard update, which leads exactly to the Ising equilibrium Gibbs measure ||]. 

Two continuous variables, x\ and A* in [—1,1], are defined on each site of a two-dimensional square lattice with 
periodic boundary conditions. Their evolution rule reads: 

^ti" = T+i- «. + 1) - 1' -i< A, < Kv 

*+i _ 2/ ( t _lA , 1 At ^ ^t /I 1^4) 



X. 



1, A* , < x\ , < 1, 



where the (time-dependent) slopes A* ■ of the Bernoulli maps are calculated according to: 

Al = tanh {J {alzl + al'^ + cr^ji, + aljl,)} , (25) 

with J a coupling constant. Discrete spin variables can be defined by: 

al=sign{xiy~xl)e{-l,l}, (26) 

which allows to retain definition Hq for the fluctuating magnetization. 

When sites are updated successively on two checkerboard lattices defined by the parity of z -I- j, it is possible to 
show that the invariant measure of this system is the same as that of the Ising model Q : 

P""^ ({o-,;j = TOij}) - exp J \^ mij{mi^i,j + rrii+i^j + mij-i + Wij+i) > ■ (27) 

Therefore, static exponents are known exactly (Onsager solution). This was checked numerically in |8|, following the 
same procedure as that used for the Miller-Huse model. Sakaguchi's model does indeed belong to the Ising universality 
class for static critical exponents. 

However, nothing is known a priori on dynamical properties, which we will investigate now along the same lines 
as before. One advantage of Sakaguchi's model is that the critical point Jc is known exactly. One source of error is 
thus eliminated. One drawback is the following: the magnetization M depends on site values at two consecutive time 
steps. The method used previously in order to prepare initial conditions with fixed small initial magnetization cannot 
work. Instead, we applied the following procedure: start from random initial conditions and evolve the system at high 
temperature (at low Jinit, the system is in the disordered phase with zero mean magnetization) until it reaches by itself 
the desired initial magnetization toq. The system is then quenched to the critical coupling. Since the magnetization 
is calculated from discrete spins, there is no dispersion on the value of initial conditions (ttiq = 2K/L'^, as before). 
Yet, this procedure is more costly numerically, which explains why, in the following, the obtained statistical accuracy 
is somewhat lower than for the previous two models studied here. In practice, we use Jinit = 0.1, before quenching to 
Jc = ln(l + V2)/2. 

The protocol used is the same as before. As expected, this model also exhibits a well-defined initial critical slip 
regime. In fact, the microscopic time imic turns out to be equal to (Fig.o). A possible explanation for this is that 
the macroscopically-correlated domains present at Ji„it survive the quench to jc- Impeccable scaling is observed for 
the initial growth of the magnetization M. Reliable measurements are possible only for L = 32 and 64: cross-over to 
the nonlinear relaxation regime occurs too soon {to ~ 60) for L = 16, and L = 128 is too costly to reach satisfactory 
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statistical quality due to the preparation of initial conditions. Ensemble average is done over 128000 realizations 
simulated during T = 128 timesteps. From linear fits over the full interval < t < 128, we obtain 9' = 0.197(7) for 
L — 32 and 0.19(1) for L = 64 (Fig.0). Within error bars, no finite-size effects are present. Our estimate is: 



-^Sakaguchi 



0.194(14), 



(28) 



in good agreement with results obtained for Model A |lq,h9|, and not consistent with the value obtained for the 
Miller-Huse model. 
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FIG. 6. Measure of the dynamic critical exponent 9' for Sakaguchi's model. The magnetization M{t), measured at the 
critical point J = Jc (cf. text), is plotted vs. time t in log- log scale. In both graphs, the solid line corresponds to a slope 
equal to 0' = 0.193 (Model A). Graph (a): the system size L = 32 is fixed, while the initial magnetization rno = IKjl? varies 
between mo = 7.8 10~^ and 9.8 10"'', for values oi K = 4,6,8, 10 (from bottom to top). The microscopic time tmic is equal to 
zero. Scaling is thus observed over more than two decades. Graph (6): the initial condition A' = 4 is fixed, for system sizes 
L = 16,32,64 (from top to bottom). The crossover time for the smallest size L = 16 is observable: to ~ 60. Finite-size effects 
are negligible for L = 32, 64. 

Again, measuring z is more difficult. The scaling of M*-^) starts late, with a fairly large microscopic time (imic ~ 40 
for L = 64, tmic ^ 70 for L — 128), estimated as before from the evolution of local exponents. As for non-zero 
initial magnetizations, ensemble average is performed over 128000 realizations, and the simulation time is T = 128, 
for system sizes up to L = 128. Exponents are thus estimated over less than one decade. However, finite-size effects 
are negligible for L > 32. We obtain: 



Csakaguchi — 0.82(2) 



(29) 



(note the large error bars). From the exact value (ijv — 1/8, this leads to: z — 2.13(5), roughly compatible with 
both Model A and the Miller-Huse model. A longer simulation time (T = 1024), even for the largest system size 
(L = 128) leads to qualitatively identical results, without any quantitative improvement: the insert of Fig. ^.a shows 
that fluctuations of the local exponent Cioc, computed over time intervals of duration tioc = 100, include values of C, 
expected for both Model A and the Miller-Huse model. 

The situation concerning the scaling behavior of the time autocorrelation function is somewhat more satisfactory, 
despite (or thanks to) finite size effects: large size behavior converges toward a value of b which is necessarily smaller 
than the Miller-Huse exponent 8 — 0.81(1). In fact, data obtained for the largest size considered [L — 128, T — 128) 
is characterized by a scaling exponent 8 = dj z — 9' ^ 0.75, compatible with Model A value 6 ~ 0.74(1) (cf. Fig. H.b). 
This finding is confirmed qualitatively by runs performed for longer simulation times [L = 128, T = 2048): the 
Miller-Huse value is not allowed by the evolution of effective behavior (exponent) with system size. Since we cannot 
cross-check with larger system size due to numerical limitations, we do not provide error bars on our estimate: 



-^Sakaguclii 



0.75. 



(30) 
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FIG. 7. Measure of the dynamic critical exponent z for Sakaguchi's model. The second moment of the magnetization M^^' (t) 
(graph (a)) and the temporal autocorrelation function A{t) (graph (6)) are plotted in log-log scale vs. time t, as measured at 
the (exact) critical point J = J^. for zero initial magnetization mo = 0, up to time T = 128. The system sizes considered are 
L = 16,32,64,128, for 128000 realizations. Graph (a): finite-size effects become negligible for L > 32. The solid and dashed 
lines respectively correspond to slopes equal to C = 0.817 (Model A) and (" = 0.839 (Miller-Huse model). In the insert, we plot 
the local exponent Cioc(^) computed over intervals iioc ~ 100, for L — 128, T = 1024, data averaged over 32000 realizations. 
Fluctuations of (ioc{t) encompass both values expected for Model A and the Miller-Huse model. Graph (6): the slope of the 
solid line is equal to —5 = —0.74 (Model A). The insert shows, for L — 128, a log-lin plot of the ratio A{t)/t~ vs. time t for 
the numerical values 5 = 0.74 (Model A, bottom curve) and 5 — 0.81 (Miller-Huse model, top curve). This suggests that the 
continuous transition of Sakaguchi's model belongs to the universality class of Model A for the exponent 5. 



Our numerical results show clearly that the critical dynamics of Sakaguchi's model belongs neither to the universality 
class of Miller-Huse model, nor to that of the sequentially updated model studied in Sec. IV A. It is moreover likely 
that Sakaguchi's model belongs to the universality class of Model A, for both statics and dynamics: the checkerboard 
update of Sakaguchi's model may thus be analogous to that of, say, Glauber dynamics of a spin system. 

Our results may be summarized as follows (see Table I for numerical values of all exponents estimated so far): 
update is a relevant parameter for dynamic critical exponents of Ising-like transitions of coupled chaotic map systems. 
Not only models with static exponents outside the Ising universality class (Miller-Huse) have also their dynamic 
exponents different from those of Model A, but models within the Ising static class (sequentially-updated Miller- 
Huse, Sakaguchi) may also have non-Ising dynamic exponents (sequentially- updated Miller-Huse). In some sense, this 
is not too surprising since update is already known to be a relevant parameter for the dynamical exponents of spin 
systems in the Ising static class. In the next section, we go a step further and investigate whether the universality 
found at the static level in S among synchronously-updated CMLs subsists at the dynamic level. 



V. NON-UNIVERSALITY WITHIN SYNCHRONOUSLY-UPDATED MODELS 



In this section, we study two variants of the Miller-Huse model: we consider first the case of a locally anisotropic 
coupling to three neighbors, and then the case where the piecewise linear local map dl3) is replaced by a smooth map. 
The static exponents /3, 7 and v of these variants were measured to be the same, within error bars, as those of the 
Miller-Huse model. Given the difficulties encountered above when trying to estimate the exponent z, in this section 
we restrict ourselves to measurements of the exponent 6' , via the early-time scaling of M{t). 

The three-neighbor variant of the Miller-Huse model was studied in [|| because of its particularly weak corrections 
to scaling. Its anisotropic evolution rule reads: 



n*+i 



.*+' 



i 



{l~ig)f{:A^.,)+g (/( 

(1 - 3g) /(a;L+ij 






.) + fi 



2i+l,j 



+ f{4^,, + l)) 



(31) 



where / is the original piecewise linear map ([l^). Each site is thus coupled to three of its nearest neighors: sites 
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belonging to even (resp. odd) columns of the lattice are coupled vertically to their northern (resp. southern) neighbor 
only. Rule ( pl|) is applied synchronously to all sites, with periodic boundary conditions. The local anisotropy 
introduced vanishes at large scales. In this case, an Ising-like transition takes place for a critical coupling gc — 
0.25118(4), estimated using Binder's method ^. 

The smooth map variant only differs from the Miller-Husc model by the choice of the local function /, which remains 
a chaotic, odd map of the [—1, 1] interval for symmetry reasons, but now reads: 



f{x) = 3x - 4x^ 



(32) 



On general grounds, this smooth function, with its expanding and contracting parts, may be considered more "generic" 
than the original piecewise linear map. An Ising-like transition also takes place in this case, for a slightly smaller 
critical coupling than the Miller-Huse model: Qc = 0.17864(4) ||g]. 

For both variants, hereafter referred to as "MH3" (three- neighbor coupling) and "C4" (smooth cubic map), the 
methodology used in Sec. Ill to determine exponent 9' for the original Miller-Huse model can be applied. Initial con- 
ditions with fixed, small magnetization mg are prepared using the same procedure, and the phenomenology observed 
is the same. In both cases, the microscopic time is visually estimated to be tmic '^ 5. After this time, a clear scaling 
behavior is observed for M{t) (Fig. ||). Sizes L = 16, 32, 64 and 128 were studied, with ensemble-averaging over 
512000 realizations for L < 64 and 128000 for L — 128. Finite-size effects are smaller than error bars for L > 32. 
Statistical errors on 6' are estimated by comparing values obtained for different initial magnetization toq = 2K/L^ 
with K = 2,4,8 and coupling strengths g within the uncertainty interval of g^- Our global (conservative) estimates 
are: 



l9(MH3) = 0.165(8) and 



'(C4) 



0.128(6) 



(33) 



Not only these values are not compatible with each other, but both are also at odds with the value found for the 
original Miller-Huse model {6' = 0.146(9)). 
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FIG. 8. Measure of the dynamic critical exponent 6' for two variants of the Miller-Huse model. The magnetization M{t), 
measured at the critical coupling gc (cf. text), is plotted vs. time t in log-log scale. The system size is L = 64 (top three curves 
on each graph), and L = 128 (bottom two curves, shifted down for clarity). The initial magnetization is mo = 2K/L^. The 
microscopic time tmic is equal to 5. Graph (a): anisotropic three-neighbor coupling with original piecewise linear local map 
( [L3| ) (model MH3). The solid lines correspond to slopes equal to 9' — 0.165. From bottom to top: L — 128, K — 8,16 and 
L = 64, K — 2,4,8. Graph (b): symmetric four-neighbor coupling with smooth cubic local map (B3) (model C4). The solid 
lines correspond to slopes equal to 0' — 0.128. From bottom to top: L — 128, K — 8,16 and L = 64, K — 3,4, 8. 

Numerical estimates of 9' obtained for all models we considered are gathered in Table II. Except in the presence of 
finite-size and finite-time effects not detectable in the experimental conditions of our work, the above results lead to 
conclude that CMLs of the Miller-Huse type exhibit different exponents for the dynamic critical properties of their 
Ising-like transitions even though they share the same static exponents and the same mode of update. 
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VI. DISCUSSION 

The first conclusion about the series of numerical experiments conducted in this work is that non-trivial early- 
time critical dynamics occurs in phase transitions of far from equilibrium, coupled map lattices. The microscopic 
time imic after which scaling sets in tends to be small for M(t), generally of the same order of magnitude as for 
equilibrium systems. On the other hand, imic is typically at least one order of magnitude larger for the evolution 
of second-order quantities such as M(^)(i) and A{t). This is also in agreement with spin systems [|l9|. One would 
however like to understand why t,„ic is always larger for second-order quantities than for M(t). Intuitively, it can 
be argued that in order for the scaling regime of M'-^^t) and A{t) to set in, the system first needs to generate 
macroscopically-correlated regions, i.e. M{t) must already be in its scaling regime. Thus t^ic{M'^'^\t),A{t)) > 
tmiciM{t)). Further, one may argue that macroscopically-correlated regions for second order quantities can only 
establish themselves when the magnetization M{t) has changed "significantly" , say by a factor 2. Using M{t) ^ t^ , one 
obtains t^ic{M^'^' , A{t)) ~ 2^'^ . This admittedly rough argument may explain semi-quantitatively why microscopic 
times observed for the Miller-Huse models with sequential and synchronous update differ by two orders of magnitude: 
tmic(Seq.)/imic(MH) - 2(i/''06)-(i/"i'*6) - 2^°, or 1024. In addition, finite-size corrections to scaling are also typically 
much larger for second-order quantities than for M{t). A direct consequence of the observed values of imic and of the 
strength of finite-size corrections is that, in practice, the exponent 9' is much easier to measure than the exponent z, 
our initial goal. Thus, although z remains largely out of reach in some cases, the critical initial slip method offers the 
advantage of accurate estimations of 9' . 

The second conclusion of our work lies in the interpretation of the various dynamic exponent values we measured, 
which are summarized in Tables I and II. 

First, and perhaps most importantly, we confirm, at the level of dynamic critical properties, the results obtained in 
Pi at the static level: update is a relevant parameter for universality classes of Ising-like transitions. Synchronously- 
updated coupled map models — such as the Miller-Huse model and its variants MII3, C4 — have critical exponents 
different from those of Model A both at the static and dynamic level. For the original Miller-Huse model, our data is 
even conclusive for both exponents 9' and z. 

Secondly, we extend to the transitions of coupled chaotic maps a result already known for equilibrium spin systems: 
the static universality class of the Ising model breaks down at the dynamic level. The sequentially-updated Miller- 
Huse model and the Sakaguchi model, which are in the static Ising class, possess different dynamic exponents (in fact 
our data is unambiguous only for 0'). The distinction between checkerboard and sequential update is known to be 
irrelevant at equilibrium for both static and dynamic critical exponents, including the exponent 9' [^Jl9|. However, 
the same distinction becomes relevant at the dynamic level for lattices of coupled chaotic maps. An additional 
interesting point is that Sakaguchi's model seems to possess the dynamic exponents 9' and z of Model A. 

Thirdly, our data gathered for the sole exponent 9' indicates the splitting of the "non-equilibrium universality 
class" for the static critical properties of synchronously-updated coupled map systems showing Ising-like transitions. 
Variants of the Miller-Huse model, all shown in B to share the same static exponents, exhibit different values of 9' . 
Again, as for spin systems, universality classes are narrower for dynamical critical exponents than for static ones. 
Excluding Sakaguchi's model, whose control parameter is too markedly different to allow a meaningful comparison, 
one may notice that numerical estimates of 9' for the four other models depend monotonously on the critical coupling 
gc (see Table II): lower values of the coupling constant correspond to a slower coarsening process. We believe this can 
be understood by considering the structure of the phase space of these systems which can be seen as a hierarchy of 
repellors. During coarsening, this hierarchy is explored, and thus its scaling properties — which probably depend on g 
— must be related to the scaling behavior which defines 9' . Of course, the above picture will need to be substantiated 
in the future, but it already suggests a link between the structure of the phase space (and hence the details of the 
dynamics) and exponent 9' . 

In agreement with our findings, which confer a sort of "maximal" non- universality to 0', the tentative picture 
sketched above, invites comments on the relative status of z and 9' . Unfortunately, our attempts at measuring z (via 
C, and 5) have largely failed, mainly because of large values of imic- However, data obtained for the variants C4 and 
MH3 of the Miller-Huse model are compatible with the idea of a greater universality for z than for 9' . In particular, 
they suggest that variant C4 might possess the same z value as the original Miller-Huse model. Thus, we would like 
to suggest that 9' may be considered as a quantity highly dependent on the details of the dynamics in non-equilibrium 
situations, while 2;, in comparison, is a more global quantity. The independence of exponents z and 9' is consistent 
with this conjecture. 

To summarize, our numerical experiments provide additional contradictory evidence to the conjecture of [0 in 
the case of continuous phase transitions of coupled chaotic map systems. In other words, the qualitative features of 
scaling are correctly predicted by renormalization-group methods (since, e.g., the scaling behavior predicted by M] 
is observed), but the same techniques fail to predict exponents quantitatively. The above remarks, though, call for 
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more detailed investigations of the respective role of the various dynamical exponents involved in order to explain 
the origin of the non-universality reported here. In particular, the possibly different status of z and Q' provides an 
interesting starting point on which to base further research, at both the numerical and theoretical levels. Further, one 
would also like to know whether the various persistence exponents defined and studied recently for spin systems |23| ] 
— and particularly the "global persistence" exponent defined in 124] — also pertain to continuous phase transitions of 
chaotic coupled map lattices, and if so, whether their numerical value also depends on the fine details of the coarsening 
process following uncorrelated initial conditions, similarly to the dynamic exponent Q' . This is the subject of ongoing 
research. 
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Model A 


Miller-Huse 


Sequential 


Sakaguchi 


Critical point 

e' 


In(l + V2)/2 
0.193(5) 


0.20534(2) 
0.146(9) 


0.11255(5) 
0.06(1) 


ln(l + V2)/2 
0.194(14) 


2 


0.808(6) 
2.165(15) 


0.839(3) 
2.07(2) 


0.83(6) 
2.12(15) 


0.82(2) 
2.13(5) 


5 = d/z- e' 


0.74(1) 


No estimate 


0.78(8) 


-0.75 



TABLE I. Update is a relevant parameter for dynamic critical exponents: a summary of numerical estimates of the dynamical 
exponents 9' , C,, z (obtained from C,) and 5, for three models with synchronous (Miller- Huse model, Sec. Ill), sequential (Sec. IV 
A), and checkerboard (Sakaguchi's model. Sec. IV B) update. The exponent values given for Model A are discussed in Sec. II. 
Number(s) in brackets correspond to the uncertainty on the last(s) digit(s), e.g. 0.193(5) means 0.193 ± 0.005, 



2D Ising Sakaguchi Sequential C4 Miller-Huse MH3 

Static universality class Ising Ising Ising Synchronous Synchronous Synchronous 

Critical point ln(l + y2)/2 ln(l + y2)/2 0.11255(5) 0.17864(4) 0.20534(2) 0.25118(4) 

9' 0.193(5) 0.194(14) 0.06(1) 0.128(6) 0.146(9) 0.165(8) 

TABLE II. The universality observed for the static properties of Ising-like transitions of lattices of coupled chaotic maps 
does not hold for their dynamic critical properties: a summary of numerical estimates of the dynamical critical exponents 9' . 
Number(s) in brackets correspond to the uncertainty on the last(s) digit(s), e.g. 0.193(5) means 0.193 ± 0.005, 
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